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Abstract. 

We report on the temporal and spectral characteristics of the early X-ray emission from the GRB 050822 as 
, observed by Swift. This burst is likely to be an XRF showing major X-ray flares in its XRT light-curve. The 

■ quality of the data allows a detailed spectral analysis of the early afterglow in the X-ray band. During the X-ray 

' flares, a positive correlation between the count rate and the spectral hardness (i.e. the higher the count rate, the 

, harder the spectrum) is clearly seen for the X-ray flares. This behaviour, similar to that seen for Gamma-ray 

' pulses, indicates that the energy peak of the spectrum is in the XRT energy band and it moves towards lower 

energies with time. We show evidence for the possible detection of the emergence of the forward-shock emission, 
^ ^. produced at a radius larger than 4 x 10^^ cm in the case of a CBM afterglow model (a formation region clearly 

I ' different from that producing the prompt emission). Finally, we show that the null detection of a jet break up to 

Pi To -I- 4 X 10^ s in the X -ray light curve of this XRF can be understood: i) if the jet seen on-axis is uniform with a 

large opening angle {6 > 20°); or ii) if the jet is a structured Gaussian-like jet with the line-of-sight outside the 
^ ' bright Gaussian core. 



Key words, gamma-ray: bursts - Gamma-rays, X-rays: individual(GRB 050822), energy peak, XRF, thermal 
component 

1. Introduction gests that XRFs and long GRBs share a similar progen- 

itor. Thus, it has been proposed that XRFs are simply 



an extension of the long-GRB population with low val- 
ues of the energy peak (Ep) of the prompt spectra (e.g. 
Sakamoto et al. 2005, Barraud et al. 2003). 



X-ray flashes (XRFs) and X-ray rich Gamma-ray bursts 
(XRR GRBs) first detected by Ginga and BeppoSAX (e.g. 
Heise et al. 2001) emit most of their prompt energy in 
X-rays (see Lamb et al. 2004). It has been shown that 

XRFs and GRBs share many observational properties, in- ^ ^^^er of theoretical models have been proposed to 

eluding: i) the temporal and spectral properties of the ^^P^^^^^ ^^^^^ ^""^ ^^^^^ ^'^t^'^^i^ P^^^^^^l 

prompt emission (e.g. Heise et al. 2001, Kippen et al. ^nces m the jet outflow (e.g. Mizuta et al. 2006) or m the 

2003, Sakamoto et al. 2005) ; ii) host galaxy properties (e.g. geometries between XRFs and GRBs. Thus, the "dirty 

Bloom et al. 2003); iii) broadband afterglows as observed f^'^^ball" invokes entrainment of baryonic material in the 

by Swift (Gehrels et al. 2004 - XRFs 050215B, Levan et al. ^^'^ j^*' ^suiting in a bulk Lorentz factor F « 300 (e.g. 

2006a; 050315, Vaughan et al. 2006; 050406, Romano et termer et al. 1999, Huang et al. 2002, Dermer & Mitman 

al. 2006; 050416A, Mangano et al. 2006; 050714B, Levan ^004). Mochkovitch et al. (2004) have alternatively pro- 

et al. 2006b; 060218, Campana et al. 2006b). The associa- PO^^^^ ^RB jets, in which the bulk Lorentz factor 

tion of XRFs with supernovae of type Ib/c (e.g. Tominga ^ > 300 and the contrast between the bulk Lorentz factors 

et al. 2004, Thomson et al. 2004, Watson et al. 2004) sug- coUidmg relativistic shells is smaU, can also produce 

XRFs. It has also been proposed that XRFs could simply 

Send offprint requests to: ogl9@star.le.ac.uk have an intrinsically wider jet opening angle in the case of 
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a uniform jet model, since the energy peak of GRB spec- 
tra is anti-correlated with the jet opening angle (Lamb et 
al. 2005; see also Li et al. 2006). 

On the other hand, other models simply invoke an ef- 
fect of the viewing angle. Indeed, Mcszaros et al. (2002) 
have stressed that X-ray photons could be produced by the 
view of the cocoon surrounding the GRB jet as it breaks 
out, instead of the narrow jet (also sec Zhang, Wooslcy 
& Heger 2004). Another interesting model based on the 
unification scheme of AGN speculates that XRFs could 
be the result of a highly coUimatcd GRB jet viewed off 
the jet axis (Yamazaki et al. 2003, Zhang et al. 2004). 

The recent Swift broadband observations of XRFs have 
shown a variety of temporal and spectral behaviour. In 
the case of the peculiar event 060218, it was established 
that the explosion was quasi-isotropic (Soderberg et al. 
2006). Similarly, Mangano et al. (2006) have shown that 
the jet opening angle of GRB 050416A could be much 
larger (6* > 20°) than those derived for GRBs (5 - 10°, 
e.g. Frail et al. 2001). On the other hand, the light curve 
of GRB 050315 (Vaughan et al. 2006) shows evidence for 
a possible jet break at 2.5 x 10^ s, implying a jet opening 
angle of 5°, consistent with the values derived for GRBs. 
Finally, XRF 050406 has been shown to be a burst possibly 
seen well off the axis of a structured jet (Schady et al. 
2006). 

These results suggest that the origin of XRFs is still 

not settled. It is therefore important to study the XRF- 
like events in detail to constrain their true origins. 

Here, wc report the case of a burst detected by the 
Swift BAT (Burst Alert Telescope; Barthelmy et al. 2005) 
on 22nd August 2005. The X-ray light curve of this event 
exhibits a stecp-to-flat-to-stcep decay, and large X-ray 
flares are superposed on the initial steep underlying de- 
cay. We show that the spectrum of one of the X-ray flares 
could present a quasi-thermal component. The paper is 
organised as follows: in Section 2, we present the char- 
acteristics of the observations and the basic steps of the 
data reduction. In Section 3, we present the temporal and 
spectral analysis of the multi- wavelength observations. We 
establish that this burst is probably an XRR GRB or an 
XRF, by using the BAT spectral results to compute the 
softness ratio (e.g. Lamb et al. 2004). In Section 4, we in- 
vestigate the physical mechanisms producing the spectral 
and temporal characteristics of the burst. 

By convention, we note hereafter the flux in the X-ray 
band is modelled as F^, ex. v~^{t — To)~", where j3 is the 
energy spectral index, a is the temporal index, and Tq is 
the BAT trigger time. We use the symbol F to refer to 
the bulk Lorentz factor. The BAT spectral slope is noted 
as /3bat- All the time intervals are hereafter referenced to 
the BAT trigger time. 



2. Observation and data reduction 

2. 1 . BA T observa tions 

The burst 050822 (Swift-BAT trigger 151486) was de- 
tected at 03:49:29 UT on 22nd August 2005 at (J2000) 
RA= 03''24™19" and Dec=-46<i01'22", with an uncer- 
tainty of 2 arc-minutes (Blustin et al. 2005). 

The BAT spectra and light curves were extracted using 
the BAT analysis software (build 2.3) as described in the 
Swift BAT Ground Analysis Software Manual (Krimm, 
Parsons & Markwardt 2004). 

2.2. XRT observations 

The X-Ray Telescope (XRT, Burrows et al. 2005) started 
to observe the burst 95 s after the trigger, following the 
sequence of readout modes: Image mode (IM) at the end 
of the slew, Windowed Timing (WT), and then Photon 
Counting (PC) modes while pointed at the target (Hill 
et al. 2004, 2005). The Low Rate Photo-Diode mode is 
no longer used, since the XRT CCD detector was dam- 
aged by a micro-meteoroid on 27th May 2005 producing 
several bad columns (Abbey et al. 2005). The XRT obser- 
vations are summarised in Table 1. Note that the X-ray 
light curve needed to be c;orrected for the loss of counts, 
because the source is located on the CCD chip close to 
the bad columns. To do that, we fitted the profile of the 
XRT point spread function (PSF, Moretti et al. 2005) to 
estimate the fraction of lost counts in the IM data and the 
first orbit of the WT and PC data (i.e. before Tq-MOOOs). 
The correction factor applied to the X-ray light curve be- 
fore To + 1000 s is / - 1.22. 

An uncatalogued X-ray source was identified at 
(J2000) RA =03'>24'^27.26" and Dec =-46^02'00.3" with 
an uncertainty of 1.4" at a 90% confidence level. This 
refined ground-calculated position was obtained after as- 
trometry corrections. To do this, we remove the first 100 s 
of each orbit where the star tracker attitude was less sta- 
ble. For GRB 050822, this leaves only 216 ks of PC data. 
The data are further filtered to remove any remaining hot 
pixels that are not filtered out by the normal pipeline 
processing, then exposure maps are made based on the 
remaining data, and all images and exposure maps are 
summed. Wc obtain all of the optical objects within 15' 
from either SDSS if available or USNO-Bl if not. In the 
case of GRB 050822 we use USNO-Bl. To find serendip- 
itous X-ray sources for matching, we run wavdetect on 
the combined XRT image, and then run xrtcentroid to 
get the best positions taking into account the instrument 
PSF and exposure maps. We do not do individual object to 
object matching, but rather we match all X-ray sources to 
all optical sources and grab all matches with a separation 
of less than 20". We look for clustering in those matches 
to find the overall mean frame shift. We find the weighted 
mean frame shift measured from all the matches and re- 
move all outliers further away then 2-cr from the mean. 
We then iterate finding the mean and removing outliers 
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for a few more iterations also requiring only one match 
per X-ray source on the third iteration. Finally, we take 
this mean shift and apply it to the GRB position. We 
calculate the statistical position errors using the empir- 
ical fits as described in Moretti et al. (2006), assuming 
that the astrometry correction removes the 3.5" system- 
atic error normally applied to XRT positions to account 
for errors in the star tracker attitude solution. We add the 
statistical error to the error from the frame shift due to 
the counting on each individual serendipitous source. We 
note that our best XRT position is 0.5" away from the 
astrometry corrected position (RA(J2000)=03''24™27.22% 
Dec(J2000)=-46<^02'00.0" with an uncertainty of 0.7" at 
90% confidence level) given for this burst by Butler (2006). 

The XRT data were processed by the Swift software 
version 2.5 ^. This software release includes new response 
files for the PC and WT modes which significantly im- 
proves the spectral response at low energy (below 0.7 keV). 
It is now possible to extend the fits down to 0.3 keV in 
both modes (Campana et al. 2006a). The residuals below 
0.6 keV are better than 10% and the flux accuracy between 
the PC and WT modes are better than 5%. A cleaned 
event list was generated using the default pipeline, which 
removes the effects of hot pixels and the bright Earth. 
From the cleaned event list, the source and background 
spectra were extracted using xselect. 

Due to pile-up in the IM data, only the WT and PC 
data were useful for spectral analysis. The PC data from 
~ 354s to ^ 414s and from ~ 617s to ~ 690s with a 
count rate above 1 count s~^ are moderately piled- up. 
The innermost four-pixel radius was excluded, and the 
source and background spectra were extracted using an 
annular region with an outer radius of 20 pixels. The same 
annular region was used to correct the pile-up effect in 
the PC data for the X-ray light curve. For the spectral 
and temporal analysis, we used the grade 0-12 events for 
the PC mode and the grade 0-2 events for the WT mode, 
giving slightly higher effective area at higher energies. The 
ancillary response files for the PC and WT modes were 
created using xrtmkarf. 

2.3. UVOT and other optical observations 

UVOT (UV-Optical Telescope; Roming et al. 2005), which 
began to observe 138 s after the trigger, detected no optical 
fading source down to a 3(7 limiting magnitude of 19.5 in 
V-band for a 278 s exposure and 19.4 in U-band for a 188 s 
exposure (Page et al. 2005). 

ROTSE-III (Rykoff et al. 2005) started to observe 
31.7 s after the trigger (i.e. during the Gamma-ray prompt 
emission phase), but no source was detected down to an 
mifiltered magnitude of 16.6 in a 84 s (at Tq + 31.7 s) co- 
added images exposure, and 17.5 in a 246 s (Tq -1-412.6 s) 
co-added images exposure. 

No redshift information is available for this event. 



^ See http;//heasarc. gsfc.nasa.gov/docs/swift/analysis/ 
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Fig. 1. 0.3-10 keV X-ray light curve of the X-ray source close to 

the position of GRB 050822. The light curve is extracted from 
~ 5 X lO'* s to ~ 5 X lO*" s using a 15 pixel-radius extraction 
region. 



3. Data analysis 

All the errors cited below are given at a 90% confidence 
level for one parameter of interest {i.e. Ax^ = 2.706). 

3.1. Spatial Analysis 

A faint X-ray source (SX) near the X-ray counterpart of 
GRB 050822 was detected at (J2000) RA=03'^24"22.37« 
and Dcc=-46''02'10.55" with a 90% error radius of 
2.3" using the same method described in Section 2.2. 
This source is 4.9" away from the X-ray counter- 
part of GRB 050822. A possible optical comitcrpart 
to SX (1.5" away from the SX position) is found in 
the catalogue USNO-Bl.O (RA(J2000) = 03^24'"22.37% 
Dec(J2000)=-46<*02'12.0"; Monet et al. 2003). 

Fig. 1 shows the X-ray light curve of this nearby source 
in the 0.2-10 keV energy band from ~ 5 x lO^s to ^ 
5 X 10® s using a 15 pixel-radius extraction region. With 
a mean count rate of ^ 10^^ counts s~^, the source does 
not contaminate the light curve of GRB 050822 before ~ 
To + 3 X 10^5 s. While counts from GRB 050822 dominated 
the field, a 20 pixel radius was used to extract the light 
curve and spectra. After ro-|-5 x 10'* s, a 15 pixel extraction 
radius was used. 

3.2. Light curve 
3.2.1. Gamma-ray band 

GRB 050822 exhibits a complex multi-peaked light curve 
(see Fig. 2), with peaks at Tq+ ~ Os, ~ 42 s, ~ 48 s and 
~ 55 — 60 s. A small peak between ~ 100 s and ~ 104 s can 
be also seen in the BAT light curve in the 15-25 keV and 
25-50 keV energy bands. A faint tail or flare from ~ 104 s 
to ~ 200 s can be also seen by eye (see the small window 
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Fig. 2. Background-subtracted BAT light curve in units of 
count s~^ (fully illuminated detectors)"^ for 4 different energy 
bands from the top to the bottom: 15-350 keV, 15-25 keV, 25- 
50 keV, 50-100 kcV and 100-350 kcV. The bin time is 1 s. The 
dashed vertical lines delimit the temporal interval where the 
BAT and XRT data overlap. In the small window using a bin 
time of 10 s, we show the late small flare between 100 s and 
104 s and a possible tail up to To -|- 200 s. 



in the top panel in Fig. 2). Above 100 keV, only weak 
emission is seen in the BAT hght curve. 

Tza and Tgo in the 15-350 keV band are 43.9 ± 0.2 s 
and 104.7 ± 0.4 s, respectively. 

3.2.2. X-ray band 

GRB 050822 shows a complex XRT Ught curve in the 0.3- 
10 kcV energy band (see the top panel in Fig. 3). The first 
1000 s of data display at least three major X-ray flares 
peaking at ~ 131s (Fi), ~ 236s (F2) and ~ 420s (F3) 
superposed on an underlying decay. 

Hereafter, we use the symbols F„ to refer to these X- 
ray flares. Note that the rise of the flare F3 is relatively 
fast with a timescale less than 30 s. After ~ 800 — 900 s, 
the light curve shows a flat-to-steep decay similar to that 
seen in other Swift bursts. 

The bottom panel of Fig. 3 shows a hardness ratio, 
defined as the ratio of the 1-10 keV band to the 0.3-1 kcV 
band, as a function of time. Some spectral hardening and 
softening are clearly seen during the rising and decaying 
parts of the X-ray flares, respectively, for the flares peak- 
ing around Tq -|- 236 s (F2) and 420 s (F3), and a clear 



GRB 050822 

Table 2. Summary of the Gamma-ray spectral fitting param- 
eters. 



Time interval 




X^dof) 


Fluence* 


since To 






(xlO"'^ erg cm-2) 


- Tso 


1.36 ±0.16 


47 (56) 


1 A+^'-i 
^■^-0.3 


T50 - 110 s 


1.63 ±0.25 


49.4 (56) 


0.9 ±0.2 


- 110 s 


1 46+°-^^ 


55.5 (56) 


^•^-0.3 


90 - 125 s 


-0.57 


53.4 (56) 


0.2 ±0.1 



*The fluences are given in the 15-150 keV energy range. 



spectral softening is seen for the decay of the X-ray flare 
peaking at Tq + 131s (Fi) (the observation began during 
this flare). 

We note that the X-ray and Gamma-ray emission 
during the temporal BAT/XRT overlap are likely to 
be produced by the same mechanism (see Fig. 3 and 
Section 3.3.2). These clues suggest that the global de- 
cay seen before 1000 s could be the result of curvature ef- 
fect emission (e.g. Kumar & Panaitescu 2000 and Dermer 
2004) associated with each X-ray flare. This is confirmed 
by Liang et al. (2006; see Figure 1 and also Table 1 in 
their paper). 

After 800-900 s, the light curve can be described by 
a broken power-law with an initial shallow slope ai = 
0.45tEj:ii followed, after a break at Utgi x 10^ s, by a 
steeper slope {a2 = 1.05 ±0.05) with a possible late X-ray 
flare peaking at ~ 1.1 x 10^ s (see Fig. 3). 

3.3. Spectroscopy 
3.3.1. Gamma-ray band 

The BAT spectra are well fit by a single power-law. All the 

spectral parameters and the fluence for different time in- 
tervals are summarised in Table 2. The use of a Band func- 
tion (Band et al. 1993) or a cutoff power-law model did 
not significantly improve the fit. The 15-150 keV fluence 
GRB 050822 is 2.3tg j x 10"^ erg cm'^ over Tgo, which 
is moderate when compared to the average BAT fluence 
of 3.1 X 10~^ erg cm^^ for GRBs from January 2005 to 
September 2006. The spectral slopes of GRB 050822 are 
relatively steep with respect to the average BAT spectral 
slope of 0.8 (e.g. O'Brien et al. 2006). This suggests 
that Ep may be below the BAT energy band, and that 
this burst may be an XRF or cLii XRR, GRB. 

Classification of a burst as an XRF or an XRR GRB 
depends on the softness ratio of the 2-30 keV fluence over 
the 30-400 keV fluence (e.g. Lamb et al. 2004). Bursts with 
SR > are classified as XRFs, bursts with -0.5 < SR < 
are classified as XRR GRBs, and those with SR < — 0.5 
are classified as normal GRBs. Because the BAT is not 
sensitive over this entire energy range, we compute SR = 

g^'^3o-4oo^kcV) integrating the best-fit spectra over 
these energy ranges. 
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Fig. 3. Background subtracted XRT light curve of GRB 050822 in the 0.3-10 keV energy band in units of flux (top panel): IM 
data (triangle); WT data (crosses); PC data (diamonds). The upper limits are given at 3 a. The times of the early X-ray flares 
arc shown in the inset. The hatched area corresponds to the temporal BAT/XRT overlap. The arrow around 1.1 x 10^ s shows a 
possible late X-ray bump (see text for more details). (Bottom panel) Hardness ratio of GRB 050822 of the 1-10 keV band over 
the 0.3-1 keV band as a function of time. The error bars are Icr statistical errors. The diamonds correspond to the PC data 
and the crosses to the WT data. The vertical dotted line represents the breaJj seen in the X-ray light curve (see Section 3.2.2). 
The dashed line is the mean value of the hardness ratio beyond To -|- 1000 s. Note that for the HR plot, we did not include the 
piled-up PC data from 354 s to 414 s and from 617 s to 690 s. 



If wc assume that the energy peak is below 2 keV, 
we find SR ^ 0.55+0 49, corresponding to an XRF. On 
the other hand, if we assume Ep = 15 keV and /?Band = 
(which is the mean value of the low energy spectral index 
of the Band function for GRBs and XRFs; Preece et al. 
2000, Kippen et al. 2003), we find SR = 0.12tJ5;?^, in 
the XRF-XRR range. A value of /^Band varying between 
— I and ^ (i.e. the range of the low-energy spectral index 
expected if the radiation is produced by the synchrotron 
mechanism; e.g. Katz 1994) would still give a value of 
SR > -0.5. So, the burst 050822 is likely to be an XRF 
or an XRR GRB. 



3.3.2. X-ray band 

The Galactic column density is N^''^ = 2.34 x 10^° cm'^ 
in the direction of this burst (Dickey & Lockman 1990). 
All the spectra were binned to contain more than 20 cts 
bin~^, and were fitted from 0.3 to 10 keV within xspec 
vl 1.3.1 (Arnaud 1996), except when the statistics were too 
low, and in these cases. Cash statistics were used (Cash 
1979). To model the absorption within xspec, we used the 
photo-electric absorption model (wabs). 



BAT/XRT analysis - To investigate whether the early 
X-ray emission is connected to the Gamma-ray emission, 
we fit the BAT and WT spectra from Ills to 125s with 
an absorbed power-law. We did not use the BAT data be- 
yond 125 s because of poor statistics, and useful spectral 
data started to be taken with the XRT only from Ills (see 
Section 2.2). A single absorbed power-law gives a good fit 
with a slope of /3 = 0.971^:^3 and an excess absorption 
value of Ann (z = 0) = 1.8 ± 0.4 x lO^^ cm'^ over the 
Galactic value {x^ /ly = 100/108) using a constant factor 
(/ = 0.95 ±0.10) to take into account the difference in cal- 
ibration between the XRT and the BAT (see Fig. 4). The 
use of an absorbed cutoff power-law or a broken power-law 
did not significantly improve the fit. This result is consis- 
tent with the hypothesis that the X-ray and Gamma-ray 
emissions are produced by the same mechanism during 
this time interval. 

XRT analysis - The fit of the WT data from Ills to 
616 s using an absorbed power-law with absorption fixed 
to the Galactic value is poor with x^/^ = 1743/264. 
Leaving the absorption component free significantly im- 
proves the fit with x^/*^ = 362/255, with excess absorp- 
tion of ANh{z = 0) = 1.4 ± 0.1 X 10^^ cm-2 over the 
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Fig. 4. Joint fit of the BAT/XRT spectra. The model is an 
absorbed power-law. 
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Fig. 5. Evolution of the peak energy during flare F3 as a func- 
tion of time, for a broken power-law spectral model. 



Galactic value. Adding a black-body component further 
improves the fit by (Ax^ = 124 for 2 DOF), and obtains an 
excess absorption column of ANh{z = 0) = 9.5^^ g x 10^° 
cm~^. A consistent value of excess absorption {ANh{z = 
0) = 1.2 ± 0.4 X 10^1 cm-2) is found when the PC data 
from ~ 800 s to ~ 4.1 x 10^ s are fitted using an absorbed 
power-law. Hereafter, we used two wabs models in xspec 
with one fixed to the Galactic value, and the other fixed 
at 1.2 X 10^^ cm~^, which is the most reliable AA'jj-value 
because the late PC data are not affected by spectral evo- 
lution (see the bottom panel in Fig. 3). 

We investigate the spectral evolution seen in the WT 
data by performing a careful time-sliced spectral anal- 
ysis of each flare. All the best-fit results and the spec- 
tral models used in each case are summarised in Table 3. 
Note that the unabsorbed fluxes given in Table 3 are cor- 
rected for the effect of the bad columns using a new tool 
XRTEXPOMAP 0.2.1 implemented in the version 2.4 of the 
XRT software. 

The spectra for flare Fi are well fit by an absorbed 
power-law with a steepening of the spectral slopes with 
time. This is in agreement with the hardness ratio (Fig. 3). 
Since we have some evidence that the X-ray and Gamma- 
ray emission is produced by the same radiation mechanism 
for a part of flare Fi, we considered fitting the spectra 
with a Band function (Band et al. 1993). However, the 
narrow 0.3-10 keV XRT energy band does not allow us to 
constrain the spectral parameters of the Band function 
with the available statistics. Instead, we approximated 
the Band function with an absorbed broken power-law 
model ^ , where /3i and /?2 are the low energy and high en- 
ergy spectral slopes, respectively. Using that model does 
not allow us to track the evolution of energy peak of the 



^ Note that within the XSPEC notations, the energy break 
(Eb) is one of the three parameters of the broken power-law 
model. We refer to Ei, in the text as the energy peak of the 
spectrum {Ep). 



spectrum {Ep) with time. These points, along with the ev- 
idence that the X-ray and Gamma-ray spectra from Ills 
to 125 s are likely to be produced by the same mechanism, 
suggest that Ep has probably already passed through the 
XRT energy band by Tq + 111 s. 

The WT spectra extracted for flare F2 (from 211 s to 
301 s post-burst) and for flare F3 (from 417 s to 616 s) 
are not well fit by an absorbed power-law, since x^/^^ = 
42.9/29 and x^/j^ = 678/135, respectively. The use of a 
broken power-law for the spectra of flare F3 allows us to 
track the decrease of the peak energy with time (fitting 
both (3i and /32, although they were tied to the same values 
for all the spectra - see Table 3 and Fig. 5). Although Ep 
was well determined during the bulk of this flare, it was 
not constrained at all for the PC spectrum from 625 s to 
789 s at the end of the flare. When this spectrum is flt by 
a single absorbed power-law, the spectral slope is still in- 
consistent with the /32-value derived from the earlier WT 
spectra. The spectral hardening seen after ~ To-|-700s (see 
Fig. 3) indicates that an extra X-ray emission component 
is probably present at this time and may account for the 
inconsistency in the /3 value (see Section 4.3). For flare 
F2, even assuming that /3i = (i.e. the mean Band func- 
tion low energy spectral slope for BATSE bursts), we ob- 
tained only upper limits for Ep for spectra WT7 and WT9, 
the spectrum WTIO being best fit by a single absorbed 
power-law with /? ~ 2. This /3 value is consistent with the 
/32-values of the broken power-law model (see Table 3). 
Although we were able to measure Ep only for spectrum 
WT8, it is not completely clear whether Ep varies during 
this flare. The results suggest that, like for flare Fi, Ep 
was below the XRT energy band for most of flare F2. 

We showed that the average spectrum of flare F3 (from 
417 s to 616 s) can alternatively be fitted using a black- 
body (BB) -I- power-law (PL) model, which was previ- 
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Table 4. Summary of the emitting radius of the black-body 
component for the WT data fit from 417 s to 501 s. 



Time interval 


Rx 


since To 


(xlO" DaoGpc 


417-431 s 




431-441 s 




441-471 s 


6.91°;? 


471-501 s 


9.81M 



-D20GPC = 20 Gpc ' where Dl is the luminosity distance 
of the source and 20 Gpc is the approximate luminosity 
distance of the mean Swift GRB redshift, using WMAP 
cosmology. 



ously used in the case of GRB 060218 (Campana et al. 

2006b), with kT = 0.185l|^ o!^^ kcV, (3 = 2.08l:[! i^, and 
X^/f = 147/133 . The BB flux corresponds to 62.5 ±3.9% 
of the total observed 0.3-10 keV flux. The use of that 
model allows us to track the decrease of the BB tempera- 
ture with time. There are also some hints that the power- 
law component steepens with time; however, the values of 
the spectral slope for each spectrum are consistent within 
the error bars, probably due to the poor statistics of the 
spectra at later times. The value of /? = 1.87lo;62 
Table 3 was obtained by tying the spectral slope to the 
same value for each spectrum. The PC spectrum from 
~ 355 s to ~ 415 s is well fitted with a single power-law 
with (3 = l-62to,ll [x^/i^ = 5/9); adding a BB compo- 
nent for this spectrum does not improve the fit signif- 
icantly (Ax^ = 1.4 for 2 dof). For the PC data from 
^ 625 s to ^ 789 s, if we fix the spectral slope of the 
power-law to /? = 1.87, then we can constrain the tem- 
perature of the black body (see Table 3). The use of a 
BBODYRAD model allows us to constrain the X-ray emit- 
ting radius of the BB component, which increases from 
Wx ~ 2.6l^;^ X 10" X (5^^) cm at the time interval 

417-431 s to ii^ ~ 9.8l x 10" x (2^^) cm at the time 

interval 471-501 s, where ]s the luminosity distance. 
We normalise to i?^ = 20 Gpc, the luminosity distance 
for the average redshift {z ~ 2.5) for Swift GRBs. At later 
time, the emission radius is no longer well constrained. 
The results are summarised in Table 4. 

The time sliced-analysis of the PC data beyond 800 s 
reveals that no significant spectral variation is seen around 
the break in the light curve at ~ 1.7 x 10'' s or the late 
X-ray bump around 1.1 x 10^ s (see Table 3). 

4. Discussion 

We established in Section 3 that the burst 050822 is an 
XRF or an XRR GRB. Its X-ray light curve shows a steep- 
to-flat-to-stecp decay. At least three X-ray flares peaking 
around To -h 131 s (Fi), 236s (F2) and 420s (F3) are super- 
posed on the initial steep decay. A strong spectral evolu- 



tion is observed during the flares. Flares F2 and F3 are best 
fit by broken power law spectral models, and we showed 
that the spectral softening during the decaying part of 
the flares is probably due to the shift of the energy peak 
of the spectrum {Ep) to lower energies. Interestingly, we 
found that the data for the X-ray flare F3 are also well flt- 
ted by a black-body plus power-law model as in the case 
of GRB 060218, the black-body component cooling down 
and expanding with time. 

We discuss the possible origin of the X-ray flares in 
the framework of the internal shock model, which is often 
invoked to interpret the such flares (e.g. Zhang et al. 2006, 
King et al. 2005). 

The X-ray light curve also shows a long smooth decay 
from ^ 1.7 X 10^ s to ~ 4 x 10^ s without any evidence 
for a jet break. We investigate whether it is consistent 
with the prediction of the current afterglow models. We 
also discuss the possible origins of the late X-ray bump 
around Tq -|- 1.1 x 10^ s and its implications. 

4.1. The origin of the X-ray flares 

The two early X-ray flares peaking around Tq -|- 236 s (F2) 
and 420 s (F3) clearly show a positive correlation between 
brightness and spectral hardness i.e. the higher the count 
rate, the harder the spectrum (see the bottom panel in 
Fig. 4), as found for instance by Ford et al. (1995) in GRB 
pulses. We also note that the temporal profiles of the X- 
ray flares are well fitted by common FRED pulse shape 
(see Liang et al. 2006). The spectral softening seen for 
the X-ray fiares can be explained by a shift of the energy 
peak to lower energy through the XRT energy band. This 
is clearly seen for the X-ray flare F3, for which the data 
are well fit by an absorbed broken power-law (see Fig. 5 
and Table 3). Indirect evidence for the shift of the peak 
energy to lower energies is also presented for the two other 
X-ray flares in Table 3. 

Internal shocks - The presence of the energy peak in 
the XRT energy band is consistent with the internal shock 
model, as shown by Zhang & Mcszaros (2002). In this 
model, the peak energy Ep of the synchrotron emission 
satisfles: 

Ep oc L^i^ ^^-1 

where L and 5t are the luminosity and the variability 
timescale, respectively. A smaller luminosity and/or a 
higher value of F and 5t produces X-ray flares rather than 
Gamma-ray peaks. Here, the 5t-values of the flares are 
larger than those of the Gamma-ray peaks. A longer du- 
ration of the X-ray flares is indeed expected at later times 
due to longer accretion episodes around the central new- 
born compact object (e.g. Perna et al. 2005 and Proga & 
Zhang 2006). It is not completely clear if the late F are 
higher or not. We could speculate that the late ejected 
shells interacting with a cleaner environment along the jet 
axis have higher F. It is nevertheless more likely that the 
main factor to produce a lower Ep is a smaller luminosity 
at later time. The shift of Ep in X-ray flares through the 
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Fig. 6. Evolution of the two parts of Eq. 1 as a function of the 
redshift (z) assuming that the Lorentz factor V has reached its 
coasting value. The intersections between the dotted and thick 
lines indicate the possible z-solutions for Eq. 1. The two parts 
of Eq. 1 are expressed in units of cm. 



for the prompt emission. It is possible that for plausible 
X-ray flare parameters, the photospheric thermal compo- 
nent may be as low as kcV or less (although more detailed 
modeling is needed, which is outside the scope of our pa- 
per). 

Assuming that the shell of matter moves rclativisti- 
cally, the variation of the emission radius in the "thin 
shell" case is given by: 

Ai?(^) ^ R^^^z)-R\{z) = ^\ti)-U V\t,)\ (1) 

where z and c are the redshift and the velocity of light, 
respectively. Here, we define f as the mean times of the 
time intervals 417-431 s and 471-501 s respectively. Fig. 6 
shows the two parts of Eq. 1 as a function of the redshift 
assuming that the Lorentz factor F has reached its coast- 
ing value. From the two plots, it seems more likely that 
the shell is mildly relativistic (F < 10 — 15). Otherwise, 
the solution of the above equation would require an un- 
reasonable high redshift (F = 100 would require a redshift 
much larger than 6). 



XRT energy band has also been reported in other Swift 
bursts (e.g. GRB 051117A, Goad et al. 2006). These au- 
thors also concluded that the X-ray flares arc produced 
by internal shocks. 

The low energy spectral slope = l.OSlg iJy) for 
the X-ray flare F3 is steep compared to the mean values 
of the (3i distribution derived from a sample of averaged 
time GRBs and XRFs (/3i ~ 0; see Preece et al. 2000 and 
Kippen et al. 2003). The steep /3i-value could suggest that 
the X-ray emission is not produced solely by synchrotron 
radiation, since the low energy spectral slope from shock 
accelerated electrons is expected to be between — i and 
i (e.g. Katz 1994, Cohen et al. 1997, Lloyd & Petrosian 
2000). 

Photospheric emission - In Section 3.3.2. we showed 
that the spectra of the X-ray flare F3 arc alternatively well 
fit by a black-body plus power-law model. According to 
the internal shock model, a quasi thermal spectrum is ex- 
pected to be produced by pair photospheric emission from 
an optically-thick shocked shell of matter becoming opti- 
cally thin at a radius i?^ . However, Comptonisation of the 
photospheric emission during the emergence of the spec- 
trum (Goodman 1986, Ryde et al. 2006, Thompson et al. 
2006) or a strong magnetic component could lead to a non- 
thermal tail in the spectrum (Thompson 1994; Meszaros & 
Rees 2000; Meszaros et al. 2002; Rees & Meszaros 2005). 
It is difficult to know whether or not the spectrum would 
peak in the XRT band, since it depends on the pair optical 
depth and the pair temperature (e.g. Zhang & Meszaros 
2002). We note that Pe'er et al. (2006) concluded that 
energy peaks below a few keV are not expected in that 
picture. However, the parameters used in that paper were 



4.2. The X-ray light curve before Tq + 800 s: the tail of 
the prompt emission 

It is likely that the X-ray light curve before Tq-I-SOOs is as- 
sociated with the tail of the prompt emission: i) the X-ray 
flares are likely to be produced by internal processes; ii) 
curvature effect emission associated with the X-ray flares 
can account for the underlying decay seen in the X-ray 
light curve before To + 800 s (see Section 3.1.2 and Liang 
et al. 2006); iii) the X-ray and Gamma-ray spectra from 
To -|- Ills to To -|- 125s are likely to be produced by the 
same physical mechanism (see Section 3.2.2). 

4.3. Constraints on the evolution of the afterglow 

The afterglow emergence - It is worth noting that a spec- 
tral hardening with time is clearly seen after ^ Tf) + 700 s 
in the WT and PC data (see the bottom panel in Fig. 3). 
Evidence that the X-ray continuum emission is sometimes 
harder during the shallow decay of the XRT light curves 
than during the initial steep decay has been found in sev- 
eral Swift bursts (e.g. O'Brien et al. 2006). The X-ray 
emission producing the initial steep decay and that pro- 
ducing the flat-to-steep decay were then interpreted as 
arising from different mechanisms (i.e. processes associ- 
ated with the prompt emission, as discussed in Section 4.2, 
and external forward shock, respectively). We argue here 
that the spectral hardening seen after ~ Tq-I- 700 s could be 
interpreted as the emergence of the forward-shock emis- 
sion. In the case of the "thin shell" CBM (Circum Burst 
medium), we could calculate a lower limit on the Lorentz 
factor from Eq. 10 in Zhang et al. (2006): 
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where tdec is the deceleration time (here t^ec 700 s). 
Eiso,52 = io^"c°rg ' ^^'^ isotropic energy of 

the burst, the efficiency for the conversion of kinetic en- 
ergy into gamma-rays and the CBM density, respectively. 
Prom the work of Sakamoto et al. (2005, 2006), it ap- 
pears that XRFs follow the Amati (Amati et al. 2002) 
relation (see also Lamb et al. 2005; Amati et al. 2007). 
Assuming that the burst follows this relation and < 15 
keV (see Section 3.2.1), the isotropic energy should be 
less than Eiso < 3 x 10^° (1 + zf erg. Since T in Eq. 
2 depends weakly on rj and n, we obtain a value of 
r > 30 X (1 + z)5/8. This would give a minimum deceler- 
ation radius of i?dec ~ 2 c U^^-^^ > 4 x 10^*^ (1 -h zfl^ 
cm. The inferred radius is much larger than the radius usu- 
ally thought of for the production of the internal shocks 
(10i3_xoi4 cm). The site of the emission after ro-F700s is 
then likely to be different to that producing the emission 
before To -|- 700 s according to the CBM afterglow model. 

Standard afterglow model - If the blast-wave evolution 
has already entered the slow cooling regime when decel- 
eration started (i.e. u > max{h'^,i'c) where Um and Uc 
are the synchrotron and cooling frequency respectively), 
then the temporal decay index (a) and the spectral slope 
(/3) after the break at ~ 1.7 x lO'^s, are predicted to be 
a = (3p — 2)/4 (a) and (3 = p/2 (b) according to the 
CBM model (e.g. Sari et al. 1998) and the wind model 
(e.g. Chevalier & Li 2000) , where p is the power law index 
of the electron distribution. We find p = 2.22 ± 0.18 using 
(b) andp = 2.07 ±0.07 using (a). These values are consis- 
tent with the commonly used values of p = 2.0 — 2.4 (e.g. 
Kirk et al. 2000, Achterberg et al. 2001). 

The shallow decay from ~ Tq -h 800 s to - 1.7 x 10" s 
can then be interpreted as a phase of energy injection in 
the blast-wave, possibly due to a longer activity phase 
of the central engine (such as the kinematic luminosity 
L oc t~'^) or a wide distribution of ejecta Lorentz fac- 
tors (Rces & Mcszaros 1998; Zhang & Meszaros 2001; 
Nousek et al. 2006; Zhang et al. 2006). In the case where 
V > max(^'i„, t'c); we find q = 0.32 ± 0.15, which is con- 
sistent with previously determined g-values (e.g. Zhang et 
al. 2006). 

4.3.1. Any evidence for a jet break ? 

The X-ray light curve of GRB 050822 after Tq + 1.7 x 
10'' s shows a monotonic, relatively smooth (except the 
late bump around 1.1 x 10^ s) and long decay up to Tq -|- 
4 X 10^ s. No indication of any jet break is seen. 

Fig. 7 shows the expected observed jet break time {tjet ) 
for different values of the observed energy peak {E°''^) as a 
function of the redshift z, using the relations from Amati 
et al. (2002; A02) and Liang & Zhang (2005; LZ05). Note 
that the LZ05 relation was originally established for op- 
tical breaks. However, if the jet models are correct, then 
the jet break time in the X-ray band should be the same. 
From the figure, it appears that whatever the values of z 
and Ep, a jet break is expected in the light curve within 
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Fig. 7. Evolution of the jet break time as a function of the 
redshift for three values of the observed peak energy -Ep''. The 
curves are deduced from the combination of the Amati et al. 
(2002) relation and Eq. 5 in Liang & Zhang (2005). The thick, 
solid line corresponds to the end of the XRT follow-up. 



the first 10 days after the burst. No such break is seen. 
A similar result was fotmd in GRB 050416A (an XRF; 
Sakamoto et al. 2006), for which the A02 and LZ05 rela- 
tions were inconsistent with the lack of a jet break up to 
To + 34.5 days. 

We discuss in the next Section whether the apparent 
absence of a jet break in the light curve can be understood 
in the framework of the current jet models. 

4.3.2. The jet models 

The model of the off-axis tmiform jet with the line-of- 
sight outside the jet edge (e.g. Yamazaki et al. 2003) can 
be ruled out. Indeed, this model predicts an initial fast 
rise when the emitting surface enters the line-of-sight of 
the observer followed by a rapid decay with a ~ p (e.g. 
Granot et al. 2002, 2005). This model is inconsistent with 
our data. 

The model of the two-component jet with the line of 
sight on or close to the less energetic wider beam is also not 
favoured. In such a model, it is expected that an afterglow 
rebrightening is seen when the fireball is decelerated so 
that the more energetic narrow component enters the field 
of view. The lack of any significant rebrightening feature 
suggests that the distinct two jet components as required 
by the model are not needed. 

GRB050416A also has a very long power-law decay 
in its X-ray light curve, with no indication of a jet break. 
Mangano et al. (2006) have modelled the X-ray light curve 
of GRB 050416A using two jet models: (1) an on-axis uni- 
form jet with a very wide opening angle (e.g. Lamb et al. 
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2005); (2) a structured Gaussian-like jet with the hne of 
sight outside the bright Gaussian core (Zhang et al. 2004). 
We can infer from their Figure 5 that either of these jet 
models could work in the case of GRB 050822. In the case 
of the on-axis uniform jet model, the lack of a jet break 
in GRB 050822 requires a large jet half-opening angle (up 
to 6* > 20°). 

4.3.3. Origin of the late X-ray bump around 

To + 1.1 X lO^s 

We next consider whether the X-ray bump around ^bump = 
To + 1.1 X 10^ s could be produced by external shocks. 
Indeed, it has been proposed that abrupt density fluctua- 
tions in the circumburst medium can produce a significant 
re-brightening in the GRB afterglows via external shocks 
(e.g. Lazzati et al. 2002). However, if the blast- wave is still 
in the relativistic regime, the flux at u > should not (or 
only very weakly) be aS'ected by circumburst density fluc- 
tuations (e.g. Nakar et al. 2003). Recent work has shown 
that if the blast-wave is still in the relativistic regime af- 
ter the interaction and v > Vc, then the decay slope a is 
expected to vary, but no re-brightening is expected to be 
seen in the X-ray light curves (Nakar & Granot 2006) . So 
the bump around 1.1 x 10^ s is unlikely to be produced by 
the result of the interaction of the blast-wave with some 
clouds of matter or density jumps. 

/^From the quality of the data around tbump, we could 
not completely rule out that the X-ray bump may be pro- 
duced by an inhomogeneity in the blast- wave or by energy 
injection when we compare the rising (0.2 < Str/thump < 
1) and decaying timescales ((5td/ibump ~ 1) of the bump 
with the limits given in Fig. 1 in loka et al. (2005). 

As an alternative, the bump around Tq+I.I x 10^ s may 
be interpreted as due to late internal shocks. Although 
this is unusual, other GRBs have exhibited some late X- 
ray flares (up to 10^ s) which were interpreted as due to 
internal shocks (e.g. GRB 050202B, Falcone et al. 2006 
and GRB 050724, Campana et al. 2006c). The quality of 
the PC data around the X-ray bump do not allow us to 
rule out this interpretation. 

5. Conclusion 

GRB 050822 is an XRF showing a complex X-ray light 
curve: i) an initial steep decay with three major X-ray 
flares; ii) a flat decay from Tq + 800 s to Tq -|- 1.7 x 10*; 
iii) a long and steeper decay up to Tq + 3 x 10^ s with a 
X-ray bump around Tq + 1.1 x 10^ s. 

We argue that the three X-ray flares observed during 
the initial steep decay are likely to be produced by internal 
processes, and that the global decay is likely to be the tail 
of the prompt emission. We showed that the energy peak 
of the spectrum for the flare peaking around Tq -|- 420 s is 
in the XRT energy band and shifts to lower energy with 
time. For the flares peaking at Tq + 131s and Tq + 236 s, 



we showed that Ep is likely to be close to or less than the 
lower end of the XRT energy band. 

Interestingly, the flare F3 is alternatively well fit by 
a black-body -|- power-law (BB-PL) model. We then pro- 
posed that the flare F3 may be produced by photospheric 
emission (involving Comptonisation) for a shell of matter 
moving at a mildly relativistic speed. 

We stress that the spectral hardening seen around 
^ To + 700 s (close to the beginning of the flat decay) can 
be interpreted as a clear indication of the emergence of 
the forward-shock emission. We showed that the emission 
after To + 700 s may then be produced in a site different 
from that producing the prompt emission, since the de- 
celeration radius should be larger than 4 x 10^^ cm in the 
case of a CBM afterglow model. 

The flat-to-steep decay can then be interpreted as be- 
ing the afterglow, the flat part corresponding to a phase 
of energy injection. The null detection of a jet break up to 
To + 3 X 10^ s in the X-ray light curve can be understood: 
i) if the jet seen on-axis is uniform with a large opening 
angle {0 > 20°); ii) if the jet is a structured Gaussian- like 
jet with the line-of-sight outside the bright Gaussian core. 
We note that the same models were also invoked in the 
case of GRB 050416A, which is an XRF (Mangano et al. 
2006) to explain the null detection of a jet break in the 
light-curve. In both scenarios, the late X-ray bump around 
To + 1.1 X 10'"' s could be produced by internal shocks, im- 
plying very late activity of the central source or it could be 
produced by inhomogeneity in the blast- wave or by energy 
injection. 
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Table 1. Log of the XRT observations for GRB 050822 following the XRT mode sequence. 
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Table 3. Summary of the X-ray spectral parameters for the best models. Note that we use two WABS models fixed to the 
Galactic value (A^f = 2.34 x 10^° cm~^) and to the value of excess absorption {ANh{z = 0) = 1.2 x 10^^ cm^^). The symbols 
F„ correspond to the three X-ray flares peaking around To + 131s (Fi), 236s (F2) and 420s (F3) as defined in Section 3.2.2. 
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^ F IS the uiiabsorbcd flux given in the 0.3-10 keV energy band in units of (xlO-® erg cm-2). 

* Due to the low statistics, we used the Churazov weighting function within xspec when fitting the spectrum (see Churazov 
et al. 1996) and a binning of 10 counts per bin. 

* The Cash statistic value and the PHA bins in parenthesis. 

° PL±BB corresponds to a power-law plus black body model. 
The spectral slope of the power-law component is fixed to (3 = 1.87. 



